Dissipation and memory effects in pure glue deconfinement 

E. S. Fraga a , T. Kodama a , G. Krein b , A. J. Mizher a , and L. F. Palhares a 

a Instituto de Fisica, Universidade Federal do Rio de Janeiro, 
Caixa Postal 68528, Rio de Janeiro, RJ 21941-972, Brazil 

b Instituto de Fisica Teorica, Universidade Estadual Paulista, 
Rua Pamplona 145, Sao Paulo, SP 01405-900, Brazil 

We investigate the effects of dissipation in the deconfining transition for a pure SU(2) 
gauge theory. Using an effective model for the order parameter, we study its Langevin 
evolution numerically, and compare results from local additive noise dynamics to those 
obtained considering an exponential non-local kernel for early times. 

Nowadays Lattice QCD provides a very reliable theoretical testing ground for the study 
of the phase structure of strongly interacting matter [ . Even results obtained within 
the simpler framework of pure gauge SU(N) can bring some insight to the analysis of 
experimental data from high-energy heavy ion collisions [ 2 j. Inspired by lattice results, 
one can build effective field theory models to be used in the study of the dynamics of 
phase conversion in the deconfining transition of non-abelian gauge theories. 

In particular, for a pure gauge SU(N) theory, the Polyakov loop provides well-defined 
order parameters [ El IH El E] > an d one can construct an effective Landau- Ginzburg field 
theory based on these quantities. The effective potential for T «Td has only one mini- 
mum, at zero, where the whole system is localized. With the increase of the temperature 
new minima appear: N minima for Z(N), the center of SU(N). At the critical temper- 
ature, Td, all the minima are degenerate, and above the new minima become the true 
vacuum states of the theory, so that the system starts to decay. 

In this paper we consider the case of pure SU(2). The transition being second-order 
[|7j, there is never a barrier to overcome. The process of phase transition will proceed 
through spinodal decomposition, which will have its exponential explosion retarded by 
the effect of dissipation in the medium, as was shown for the case of the chiral transition 
in Ref. Using the effective model proposed in Ref. [ Ej, we study the Langevin 

evolution of the order parameter, comparing results, for early times, from local additive 
noise dynamics to those obtained considering non-markovian effects brought about by an 
exponential non-local kernel. The choice of the kernel is such that one can trade off the 
memory integral for higher-order time derivatives in the differential equation for the order 
parameter. The simple form of the kernel does not restrict our semi-quantitative analysis 
since the only physically significant parameter in the memory kernel is its width, so that 
different choices for its functional form should not modify the output appreciably. 
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The effective theory we adopt [ E] is based on a mean field treatment in which the 
Polyakov loops are constant throughout the space and the free energy is a function of its 
eigenvalues, Pjk = exp(i9j) 8^. To a purely perturbative computation of the free energy of 
gluons in SU(N) one adds, phenomenologically, a mass scale M = M{T d ) to the dispersion 
relation, i.e. tOk — Vk 2 + M 2 . is the critical temperature for deconfinement that can 
be read off from lattice data. For SU(2), = 302 MeV [E]. Parametrizing the Polyakov 
loop as diag[exp(i(p, —i<j))] and defining a more convenient variable ip = 1 — 7T0/2, which 
plays the role of the order parameter for deconfinement, one obtains: 
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so that (tp) = in the confined phase and (?/>) — > tpo in the deconfined phase. One can 
easily connect the behavior of ijj to that of the trace of the Polyakov loop, also used as an 
order parameter viaTrL = 2cos[7r(l — ip)/2]. The value of M can be determined from 
the deconfinement temperature through the relation = (3/2) 1 ^ 2 M/n w 0.38985M, and 
the minima of the potential occur at i/)q — ± (1 — Tj/T 2 ) 1 ^ 2 . 
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Figure 1. Volume average of the SU (2) or- 
der parameter normalized by the positive 
minimum of the bare effective potential. 



Figure 2. Comparison of the early time 
(exact) evolution of the SU(2) order pa- 
rameter for different values of r. 



Let us now consider the nonequilibrium evolution of the order parameter for the spon- 
taneous breakdown of Z(2). As costumary, we assume the system to be characterized by 
a coarse-grained free energy in the Landau- Ginzburg fashion 



d 3 x 



a(T) 



(v^) 2 +v eff m,T) 



(2) 



where V e ff([ip],T) is the effective potential obtained previously, cr(T) = 7i 2 T/g 2 plays the 
role of a surface tension and g = g(T) is the gauge coupling. 
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A markovian description of the time evolution of the order parameter and its approach 
to equilibrium can be implemented through the traditional additive-noise Langevin equa- 
tion 



o{T) 



+nT)^t+v; ff m,T)=t(x,t) 



(3) 



The function £ is a stochastic noise assumed gaussian and white so that (£(x,t)) = 
and (£(x,t)£(x',t')) = 2T5(x — x')5(t — t'). Y is a kinetic coefficient that we fix in the 
following way. Performing pure-gauge euclidean lattice Monte Carlo simulations in the 
line discussed in Ref. [Ej, spinodal decomposition is obtained through local heat-bath 
updates of gauge field configurations at (3 = 4/g 2 = 3 (corresponding to T = 6.6 T<f)> after 
thermalizing the lattice at (3 = = 2. The critical value for deconfinement is found 
to be (3d ~ 2.3. T is then extracted by comparing the short-time exponential growth 
of the two-point correlation function predicted by the simulations [ IIP] to the Langevin 
description assuming, of course, that both dynamics are the same (see, also, the extensive 
studies of Berg et al. [ITT]). Assuming that typical thermalization times are of the order 
of a few fm/c, one can relate Monte Carlo time and real time and obtain V ~ 10 3 fm -2 . 

To incorporate memory effects in our description, we consider the generalized time- 
dependent Landau-Ginzburg evolution equation [ [T2J US] 



#V dip 
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For analytical purposes, it is convenient to choose for the kernel E(£ — t') = e~( t ~ t '^ T /V, 
which is a representation of the Dirac delta function in the limit r — > 0, in which the 
description becomes markovian. Different functional forms for a localized kernel that 
yields the markovian result in the appropriate limit are essentially equivalent, the only 
relevant scale being the width r -1 . For the particular choice above, one can dispose of 
the memory integral in favor of an additional third-order derivative in time, obtaining the 
following evolution equation for the order parameter: 



+ [* + tT(T)} ^ + - a(T)V 2 ^ + T) = t) 



(5) 



which reproduces (|3J) in the limit of vanishing r. 

We solve Eq. (jSJ) numerically on a lattice, average over several realizations with random 
initial configurations around ip ~ 0, and study the evolution of the volume average (tp) 
normalized by ipo. In Fig. 1 we plot (ip), normalized by the positive minimum of the 
bare effective potential, for three situations: no dissipation and no noise (dotted curve), 
markovian dissipative evolution with no noise (dashed curve) and with white noise (solid 
curve). When considering noise, we have added the appropriate counterterms to make 
the equilibrium solution independent of the lattice spacing [I14j. From this comparison, 
it is clear that dissipation and noise delay appreciably the thermalization time scale for 
the deconfining transition as was discussed in [ITB]. 

In Fig. 2, we show results for the early-time (exact) solution of Eq. (JHJ for very 
conservative values of the width r _1 . Here, we have linearized the effective equation 
of motion and solved it for the homogeneous average of ip at early times, a process 
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which allows for treating noise exactly [ ll|. From Fig. 2, it is clear that even for 
small values of the reaction time the delay for the onset of the deconfming transition 
is quite significant. Although our analysis is restricted to the short-time behavior (up 
to ~ 1.5 fm/c) of the order parameter in a pure gauge theory, it indicates that one 
should include non-markovian effects in more realistic approaches to the phenomenology 
of high-energy heavy ion collisions. For instance, large modifications brought about by 
memory kernels in a nonequilibrium calculation of dilepton production seem to play an 
important role in understanding the data [HZj. Strictly speaking, Eq. (jSJ) does not have 
a parabolic structure and may exhibit causality problems similar to the ones found in the 
usual diffusion equation [ EI] • On the other hand, the inclusion of retardation effects as 
discussed above was shown to be crucial for satisfying causality constraints in relativistic 
dissipative hydrodynamics, and may provide an alternative to a second-order formulation 
of thermodynamics [I19j. 
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